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Abstract 

The sextic oscillator is proposed as a two-parameter solvable 7-independent potential in the Bohr 
Hamiltonian. It is shown that closed analytical expressions can be derived for the energies and 
wavefunctions of the first few levels and for the strength of electric quadrupole transitions between 
them. Depending on the parameters this potential has a minimum at /3 = or at /5 > 0, and 
might also have a local maximum before reaching its minimum. A comparison with the spectral 
properties of the infinite square well and the f3^ potential is presented, together with a brief analysis 
of the experimental spectrum and E2 transitions of the ^^^Ba nucleus. 
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I. INTRODUCTION 



In the last few years there has been considerable interest in looking for analytic solutions 
of the Bohr Hamiltonian which describes the collective motion in nuclei in terms of shape 
variables (/?, 7) Q]. This was initiated with the introduction of the interacting boson model 
(IBM) j2|. This model present four different dynamical symmetries each one associated to 
a well defined nuclear shape. The model Hamiltonian provides a natural way of going from 
a phase to another one by changing systematically few parameters and, consequently allows 
to study shape phase transitions. Thus, the IBM phase diagram has been analyzed from 
different points of view 

mm. 

Specially important are the critical points since in 
these situations structural changes occur rapidly and is difficult design appropriate theoret- 
ical models. Recently, lachello has proposed a new dynamical symmetry for describing the 
critical point at the transition from spherical to deformed 7- unstable shapes Q]. This sym- 
metry has been called E(5) and is expected to occur in nuclei in which the V{j3, 7) potential 
depends only on the (3 variable, and it has a relatively flat shape in the (3 variable. Around 
this critical point of the phase transition the potential appearing in the Bohr Hamiltonian 
can be approximated by an infinite square well in the f3 variable. In this case the Bohr 
Hamiltonian can be solved exactly in terms of Bessel functions, and various quantitative 
predictions can be obtained for various spectroscopic properties (ratios of the excitation 
energies and B{E2) transitions), on the basis of which one can search for candidates for the 
5) symmetry among nuclei. Lately, other two dynamical symmetries X(5) and Y(5) 
I to describe the critical point at the transition from spherical to axially deformed shapes 
and from axially deformed to triaxial shapes, respectively, have been proposed. 

The introduction of the E(5) symmetry renewed the interest in studying exactly solvable 
V{(3,'y) potentials. Most efforts have concentrated on 7-independent potentials, for which 
the most well-known example is the harmonic oscillator in the five-dimensional space ^| 
and its extension which contains a term proportional to jl^. More recently two more 
exactly solvable 7-unstable potentials have been discussed: the Coulomb and the Kratzer 
potentials j3| • Here again the latter potential is an extension of the former in the sense that 
it contains a term proportional to which formally changes the r variable, the analogue 



of the I angular momentum of radial potentials in three spatial dimensions. T 



lis formal 



changing of r introduces a minimum of the potentials at /5 > in both cases 



laUJ. The 
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bound solutions of the harmonic oscillator and the Coulomb potentials (and their extensions) 
are given in terms of generalized Laguerre polynomials Q|. 

Similarly to the three-dimensional case, these examples, together with the infinite square 
well, practically exhaust those exactly solvable potentials that contain a term. Here we 
propose another potential which has this property, and although it is not exactly solvable 
in the classical sense, it has a number of features that make it an ideal potential to be 
used in the Bohr Hamiltonian. This is the sextic oscillator, which belongs to the class 
of ,uasi-e.actlv .o.vab.e (QES) potentials Q. These potentials have the property that 
their solutions can be obtained in closed form for a number of energy eigenvalues, i.e. for 
the lowest few values of the n principal quantum number, which are also the lowest in 
energy. This is clearly sufficient for potentials appearing in the Bohr Hamiltonian: it is 
rarely necessary to consider more than a few levels with the same angular momentum, and 
these can be obtained from the lowest few solutions of the sextic oscillator. Furthermore, 
the sextic oscillator has a more flexible shape than other solvable potentials, as depending 
on its parameters, it can have a minimum at /5 = or at /5 = /Smm > 0, and in addition, it 
can also have a local maximum at /Smax < Pmin- 

The paper is structured as follows. In Sect. 2 the Bohr Hamiltonian is revised together 
with its solutions for 7-independent potentials. In Sect. 3 the lowest energy solutions of the 
Bohr Hamiltonian for the sextic oscillator potential are worked out. Section 4 is devoted to 
show the simple use of the sextic oscillator as a flexible 7-independent potential. Finally, in 
Sect. 5 we present preliminary applications and discuss future extensions. 



II. THE BOHR HAMILTONIAN FOR 7-INDEPENDENT POTENTIALS 

Let us first consider the Bohr Hamiltonian describing the collective motion of a deformed 
nucleus in the five-dimensional space determined by the 9i Euler angles (z = 1, 2, 3) and 
the intrinsic (3 and 7 variables is Q| 

In what follows we assume that the potential in ((H) depends only on /3, i.e. y(/5,7) = U{[3). 
For these 7-independent potentials the wavefunctions can be separated into two parts 

^(/?,7,^^) = /(/?)$(7,^.) , (2) 
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which satisfy the following differential equations: 



sin 37 97' ^(97 ' 4 ^ sin^(7 — |7rA;) 



A = T(r + 3), r = 0, 1, 2, ... (3) 



Here we have introduced e = and = ^U{P). Note that the r values determine 
the allowed angular momenta J too [l^. By setting 0(/3) = we obtain an equation 

which has the form of a radial Schrodinger equation 

dV /(r + l)(r + 2) 



Note that this is different from the Eq. (6) in in that it contains no linear derivative term 
due to the different definition of This choice also implies that the factor corresponding 
to the f3 volume element in the integration of functions of the type /(/?) in (0)) has been 
transferred to the solutions of the type (f){(3) in (0). Thus, in the integration of these no 
factor arising from the volume element has to be included. The complete solution of the 
problem implies the solution of Eq. (jH)) too, this was solved in Ref. [lfi|. 

III. THE SEXTIC OSCILLATOR 

The sextic oscillator with a centrifugal barrier is defined as 

if = + ~ ^^^^i^' ~ ^^^^ + (b' - Aais + ^ + M)) + 2abx' + a'x' , (6) 

where x G [0, oo) and M is a non-negative integer. This potential is quasi-exactly solvable, 
which means that for any non-negative integer value of M, M + 1 of its solutions can be 
obtained in an algebraic way. The (unnormalized) solutions are written in the form 

0„(x) = P„(x2)(x2)-1 exp l^—x' - ^x'^ n = 0, 1, 2, . . . (7) 

where Pn is a polynomial of order n. Obviously, normalizability requires a > 0, while a = 
reduces the problem to the exactly solvable harmonic oscillator. 
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The simplest solutions are obtained for M = and M = 1 Q|. For M = only 
found-state) solution appears at = 4- 

wavefunction being 



one nodeless (i.e. ground-state) solution appears at E^^^ = Abs, with the corresponding 



0r=°H^)~(^^r^exp(-^x^-^a:^j . (8) 

For M = 1 two solutions appear, one nodeless, and another with one node for x > 0. 
These correspond to the ground-state and the first excited state, respectively, at energies 
eJ^^=i) = 46s + A_(s) and e[^=^^ = Abs + A+(s), where 

A±(s) = 26±2(62 + 8as)^/2 (g) 

are the roots of the equation — 46A — 32as = 0. The corresponding wavefunctions are 

<^r-)(a:) ~ (l - {x^r^ exp (-^x^ - ^x^) , (10) 

and the A = A_(s) and A = A+(s) choice has to be made for n = and n = 1, respectively 
jisl ]. (Note that a > and s > imply A_(s) < 0, so the polynomial part of pO|) is nodeless.) 
It has to be mentioned that the solutions for M = and M = 1 belong to different sextic 
potentials if s is the same, as the coefficient of the quadratic term is different then. We shall 
see, however, that with appropriate combinations of s and M it is possible to solve sextic 
potentials that differ only in the strength of the centrifugal term. 

The normalization of the wavefunctions can also be given in closed form. For this one 
has to evaluate integrals of the type 

/(^) = r x^exv(--x^-bx' 
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^ir(:i±i),.r^.(l±i,l;g, 

where Dp{z) is the parabolic cylinder function and U{a,(3;z) is one of the forms of the 
confluent hypergeometric function Q|. 

Larger values of M can also be considered (e.g. for M = 2 three different solutions are 
obtained for the three roots of a cubic algebraic equation for A) , but M = and M = 1 are 
sufficient for our purposes in this paper. A complete study including more solutions, explicit 
closed forms for the normalization factors and applications to actual nuclei is underway. 
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IV. APPLICATION AS A 7-INDEPENDENT POTENTIAL 



In order to cast (jH)) in a form similar to ^ we have to write x = (3 and = | + | 
(remember that r > 0). In order to keep the quadratic term at a constant value (once the 
a, h parameters are fixed) we also have to prescribe 

s + M + ]^ = ]^(^ + 2M = c = const. (13) 

With this the sextic oscillator Hamiltonian can be brought to the form of (@)) with u{j3) 
being 

u^[p) = (52 _ 4ac'^)/52 + 2a6/3^ + a^(3^ + , (14) 

where the index vr = ± is included to distinguish the potential for even/odd r's which is 
slightly different as explained below. In Eq. ()14|) are the constants obtained in Eq. ()13|) 
for even/odd values of r and we have introduced a constant Mq for convenience as it will be 
discussed below. 

Equation (fT^ implies that increasing/decreasing M in one unit has to come with de- 
creasing/increasing T in two units. Thus, once the values of the (a, 6) parameters are 
fixed, the sequence of (M, r) values {K, 0), {K — 1, 2), {K — 2, 4), . . . correspond to solu- 
tions of Eq. (fTl|l with c+ = | + -ft'. In the same way, the sequence of (M, r) values 
{K, 1), (ir - 1, 3), (fC - 2, 5), . . . correspond to solutions of Eq. (HH) with c' = 1 + K. 
Consequently, the potential for r-even and r-odd states is slightly different due to the fact 
that the coupling coefficient — Aac^ of the quadratic term is different in the two cases due 
to the choices for c"*" and c~ that are necessary to separate the (r + l)(r + 2)/5~^ term in a 
uniform way. This situation can be handled using different strategies. One possibility is set- 
ting Uq=Uq=Q and considering 6^ > iQa, which minimizes the deviation of the quadratic 
terms compared to the quartic and sextic terms. Another possibility is introducing a relative 
energy shift between the r-even and r-odd potentials by setting Uq and Uq such that the 
potential minima are at the same energy. We shall discuss this possibility after analyzing 
qualitatively the spectrum and the potential shapes. 

Let us analyze the spectrum obtained in a simple case. Taking M = 1 we obtain, as 
explained in the preceding section, two solutions for r = (one with no nodes, n = 0, and 
the other one with one node, n = 1). In the notation introduced in Ref. |j8| the label ^ 
is our n + 1. Thus the two solutions of Eq. (fTUI) with s = 5/4 are 0i*^~^''(/3) with n = 



6 



and 1 and correspond to (p^^r = 0i,o and 0^,^ = 02,o respectively in Ref. |8| notation. 
Inspecting the energy eigenvalues, the corresponding energies are Ei{M = l,r = 0) = 
Ei,o = 7b- 2(62 ^ i0a)i/2 + ^+ and E2(M = 1, r = 0) = ^2,0 = 76 + 2(6^ + lOa)^/^ ^ 
respectively. The same potential is obtained by taking M = and r = 2. In this case we 
have a single solution, Eq. (jSJ with s = 9/4, 4'q^^~^\(3) with no nodes and with enerry 
Ei(M = 0,r = 2) = 96 + ■Uj]'. This corresponds to (f)^^r = 0i,2 in the notation of Ref. jsf. 
A similar analysis can be performed for the solutions with odd r-values. For M = 1 there 
are two solutions with r = 1, Eq. (fTUj) with s = 7/4, which correspond to = and 
0C,T- = 02,1 for n = and n = 1 respectively. The corresponding energy eigenvalues are 
Ei^i = 96-2(62 + 14a)^/2 + Mo and i?2,i = 96 + 2(6^ + Ma)^/^ + Mo . Again the same potential 
is obtained for M = and r = 3. In this case there is a single solution with no nodes, Eq. 
(jHI) with s = 11/4, which corresponds to 01,3 and has an energy Ei^s = 116 + Mq . In Fig. ^ 
a schematic spectrum is shown with indication of the relevant quantum numbers. In Fig. |21 
the corresponding wavefunctions with the notation 0^,^ are presented. 

Now we analyze the different potential shapes that can be produced by different election 
of parameters in Eq. ()14|1 . From (fT^ we find that the shape of the potential u^{l3) depends 
on the sign of 6^ — Aac^ and 6, which set the coefficients of the quadratic and quartic terms. 
(The coefficient of the leading sextic term is always positive.) When 6^ > Aac^ and 6 > hold 
(i.e. for 6 > 2(ac'^)^/2), the potential has a minimum at /5 = and it increases monotonously 
with 13. When 6^ < Aad" , irrespective of the sign of 6 (i.e. for —2{ac^Y^'^ < b < 2{ac^y^'^), a 
minimum appears for /? > 0, while for 6^ > Aac^ and 6 < (i.e. for 6 < — 2(ac'^)^/2), first a 
maximum appears and then a minimum as (3 increases. In all three cases the exact location 
of the extremal point (s) can be obtained from the real and positive solutions of 

(/5o^)' = ^[-26±(62 + 12ac-)^/2] . (15) 

Due to the relatively small difference in c"*" and c~, the r-even and r-odd potentials have 
the same types of extrema at about the same (3, except for some peculiar combinations of 
a and 6. Assuming that there are no complications of this kind, we can now return to the 
question of renormalizing the minima of the r-even and r-odd potentials. For 6 > 2(ac^y^'^, 
TT = +, — the minima of the two potentials will be Uq and at /3 = 0, so they coincide if 
Uq = Uq holds. For 6 < 2(ac'^)^/2 we can equate the minima of w^i^P) and u~{[3) if we set 
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■Uq = and 

uo = {b' - llamf - {b' - Ua){f3,f + 2a6[(/5j-)^ - {(3^)'] + a%(3^r ' (/^o")'] > (16) 

where the are obtained from with the choice of the "+" sign. With this the two 
potentials have their minima at the same energy, but they take on different values at the 
origin. Illustrations of the possible potential shapes are displayed in Fig. El Obviously, Uq 
in also has to be added to the energies of the r-odd potential. Fig. HI show the relative 
position of the energy levels E^^r as either a or 6 is varied and the other parameter is kept 
at a fixed value. 

The electric quadrupole transition rates can also be determined analytically by calculating 



the matrix elements of the transition operator |8L 111] 

T(^'^ = ta,, = tf3 [Z^S^OST + 2-'/\D^^}, + D^L,) sin 7] • (17) 

The radial integrals that appear in the P variable in the matrix elements of T^^"^^ can again 
be determined using (|TT| . In order to obtain the total matrix elements, one has to calculate 
also the components depending on 7 and the Euler angles 6'j. This can be done following 
the techniques described in Ref. These parts introduce certain selection rules not only 
for the angular momenta, but also for r. 

V. DISCUSSION 

In order to compare the main characteristics of the sextic oscillator as a 7-unstable 
potential with those of other potentials of this kind, we present calculations for a particular 
value of the parameters, a = 40000 and b = 200. These numbers were chosen such that the 
resulting e nerg y spectrum approximates that of the ^^^Ba nucleus, the first candidate for E(5) 



symmetry 17[. We stress that our aim is not to reproduce the experimental data, rather 
to get a qualitative picture about the general performance of the model. The potentials 
M^(/3) are displayed in the middle panel of Fig. El while the energy eigenvalues are shown 
in Fig. El together with the corresponding experimental energy levels. Fig. El also shows 
the calculated and the experimental B{E2) values for transitions between the energy levels. 
Note that electric quadrupole transitions which change r with more than one unit are zero if 
we use the transition operator (fTTI) . but finite _B(E2) strengths can be obtained if we apply 
terms of the next order (see e.g. jlS^). 
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In Table m we summarize the ratio of the most important energy eigenvalues and those of 
the most characteristic B(E2) transition rates obtained from the sextic oscillator with pa- 

n 

rameters a = 40000, b = 200, the infinite square well potential [8|] and the numerically solved 



potential 



19(1 together with the corresponding experimental values for ^'^^Ba, whenever 



available. It is seen that the energy ratios corresponding to the E(5) symmetry systemat- 
ically fall between the values of the potential and the sextic oscillator. The situation 
is less obvious for the ratio of the B{E2) values: here the sextic oscillator and the infinite 
square well seem to yield similar ratios, while the numbers obtained from the potential are 
systematically higher. This might be due to the fact that the sextic oscillator potential goes 
to infinity steeper than the potential, so the asymptotic behaviour of its wavefunctions 
can be closer to that of the wavefunctions of the infinite square well. Comparing the results 
with the experimental data for ^^^Ba we can conclude that, at least in this case, the sextic 
oscillator allows a better approximation than the other essentially parameter-free potentials. 
We expect that this conclusion will be general due to the flexible nature of the sextic po- 
tential whose shape is governed by two parameters. In fact this potential can be used not 
just at the critical point but it can be useful to model the full shape phase transition from 
spherical to deformed 7— unstable nuclei by changing the parameters a and b. 

Before closing, we mention some aspects of the sextic oscillator that might give further 
help in the analysis of nuclei near critical points. First we note that with M = 2 in © the 
analysis can be extended to further states, such as 0i,4, 0i,5, 02,2, 4>2,3, 4>3,o 4>3,i- Second, 
the potential shape which contains both a local maximum and a minimum at /5 > might 
be useful in the description of nuclei with the so-called X(5) symmetry, which is thought to 
occur in the shape phase transition between the spherical and the axially deformed domain 
j^. Third, there are further quasi-exactly solvable potentials both with confining and non- 
confining nature which can also be considered in the Bohr Hamiltonian. 
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TABLE I: Ratios of some energy eigenvalues and electric quadrupole transition strengths from 
the sextic oscillator with a = 40000, b = 200, the infinite square well ^] and the potential 



together with the experimentally observed quantities for ^^^Ba. 








^(6^3) 

mil) 


S(i?2;4+2^2+i) 
B{E2;2+,^0+o) 


BiE2;2+,^0+o) 


B(£;2;0+3^2+2) 
BiE2;2+,^0+o) 


sextic osc. 


2.39 


3.68 


3.70 


1.70 


1.03 


2.12 


E(5) 


2.20 


3.03 


3.59 


1.68 


0.86 


2.21 




2.09 


2.39 


3.27 


1.82 


1.41 


2.52 


i34Ba (exp.) 


2.31 


3.57 


3.65 


1.56(18) 


0.42(12) 
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FIG. 1: Schematic typical spectrum for the sextic oscillator with indication of the relevant quantum 
numbers. 
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FIG. 3: Potentials (full line) and u-{l3) (broken line) for a = 40000 and h = 1000 (left 

panel), h = 200 (middle panel) and b = —1000 (right panel). The lowest energy level appears in 
these potentials at o = 4633.57, 73.35 and —9366.43, respectively. 
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FIG. 4: Excitation energies E^^ = E^^j- — Ei^q with a = 40000 fixed as a function of b (left panel) 
and with b = 200 fixed as a function of a (right panel) . 
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